Analysis of the Neutralino System in 

yo", Three— Body Leptonic Decays of Neutralinos 

o 
o 

(N 

S i S.Y. Choi 1 , B.C. Chung 2 , J. Kalinowski 3 , Y.G. Kim 4 and K. Rolbiecki 3 

Department of Physics, Chonbuk National University, Jeonju 561-756, Korea 
2 Department of Physics, KAIST, Daejon 305-017, Korea 
3 Institute of Theor. Phys., Warsaw University, Hoza 69, 00681 Warsaw, Poland 
. 4 Department of Physics, Korea University, Seoul 136-701, Korea 

<N 

o 
in 
o 



X 

5-H 



Abstract 



Neutralinos x° in supersymmetric theories, the spin-1/2 Majorana-type superpart- 
ners of the U(l) and SU(2) neutral electroweak gauge bosons and SU(2) neutral 
Higgs bosons, are expected to be among light supersymmetric particles so that 
they can be produced copiously via direct pair production and/or from cascade de- 
cays of other sparticles such as sleptons at the planned Large Hadron Collider and 
the prospective International Linear Collider. Considering the prospects of having 
both highly polarized neutralinos and possibility of reconstructing their decay rest 
frames, we provide a systematic investigation of the three-body leptonic decays of 
the neutralinos in the minimal supersymmetric standard model and demonstrate al- 
ternative ways for probing the Majorana nature of the neutralinos and CP violation 
in the neutralino system. 



1 Introduction 



The search for supersymmetry (SUSY) is one of the main goals at present and future 
colliders since SUSY is generally accepted as one of the most promising concepts for 
physics beyond the standard model (SM) pp. All SUSY theories contain neutralinos, the 
spin-1/2 Majorana superpartners of neutral gauge bosons and Higgs bosons, that are 
expected to be among the light supersymmetric particles that can be produced copiously 
at future high energy colliders. 

Once neutralino candidates have been detected at high energy colliders such as the 
Large Hadron Collider (LHC) |2j and the International Linear Collider (ILC) [3], it is of 
great importance to verify that the observed states are indeed the spin-1/2 superpartners 
of the neutral SM gauge and Higgs bosons. For that purpose, it will be crucial to measure 
their quantum numbers and to confirm that they are indeed Majorana fermions jU El El IZj- 
Moreover, masses, mixings, couplings and CP violating phases must be measured in a 
model-independent way jSJ El UD| to reconstruct the fundamental SUSY parameters and 
to verify the SUSY relations at the electroweak scale, leading to a reliable extrapolation 
to the grand unification scale or the Planck scale 

In this report we focus on probing the Majorana nature and CP properties of neu- 
tralinos in the minimal supersymmetric standard model (MSSM) through the charge (C) 
self- conjugate three-body decays of polarized neutralinos into the lightest neutralino Xi 
and a lepton pair 

t-x\ti- (1) 

with i = e or fi whose four-momenta can be measured with great precision. In particular, 
the decays of the second lightest neutralino x° wm be studied in more detail since it is 
expected in most supersymmetric scenarios [T2] that the three-body decay mode has a 
significant branching fraction only for the second lightest neutralino, while the heavier 
neutralinos decay mainly through two-body decays. 

Thorough analysis of the CP properties and Majorana nature of neutralinos produced 
in pairs in e + e~ annihilation has been performed in Refs. [3] and jH] [CP asymmetries 
in neutralino production with two-body decays have been investigated in Ref. [H3]. The 
spin-1/2 neutralinos xi are produced polarized with the degree of polarization depending 
on their production mechanism and polarization of the colliding beams. Because in general 
the momenta of final particles do depend on the neutralino polarization, full account of 
spin correlations between production and decay processes is necessary E] • 

Since the verification of Majorana character of neutralinos and their CP properties are 
of fundamental importance, they have to be scrutinized in all possible ways. Here we con- 
sider neutralinos Xo, which themselves are decay products of scalar particles, i. e. sleptons. 
In such a case their subsequent decay can be analyzed independently of the production 
mechanism making the theoretical treatment greatly simplified and the interpretation of 
various observables more transparent. 
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As noted explicitly in a recent work [Hj, neutralinos x\ produced in e\ decays are 
100% polarized, having negative helicity in — > e~x% and positive helicity in e\ — > e + x\- 
Furthermore, it is possible to reconstruct the rest frame of the neutralino X2 111 a ^ ew 
specific cascade processes, for example, in the process 

e + e - ^ g+g- ^ e + x?e~X2 (2) 

followed by the three-body decay X2 ~^ as shown in Refs. [TH] and [IB]. It is 

the purpose of this paper to show that such a perfect neutralino polarization combined 
with the study of angular correlations in the neutralino rest frame can provide us with 
alternative ways for probing the Majorana nature of the neutralinos* and CP violation in 
the neutralino system. 

Keeping in mind the above aspects, we provide in the present work a systematic 
analysis of the neutralino decay x\ ~~ * Xi^ + ^~ m ^ s decay rest frame for extracting all the 
physical implications due to the Majorana nature as well as CP violation |q1I§1 Hoi 1181 IT§] 
in the neutralino system of the MSSM assuming 100% neutralino polarization. Through 
the present work it is assumed that the neutralino masses have already been measured with 
great precision j20j- On the other hand, the efficiency of reconstructing the X2 polarization 
as well as its rest frame depends not only on the cascade processes under consideration 
and on the values of relevant SUSY parameters, but also on details of experimental setup. 
Experimental simulations with realistic reconstruction efficiencies and with background 
processes included, however, are beyond the scope of the semi-theoretical treatment in the 
present paper. Nevertheless, we hope that our findings are interesting enough to motivate 
realistic simulations. 

The layout of the paper is as follows. In Sect. El we briefly recall the mixing formalism 
for the neutral gauginos and higgsinos in the CP non-invariant theories with complex 
phases. In Section |3] the leptonic decays x® ~ * °f a polarized neutralino x°i 

(i = 2, 3, 4) are described in terms of quartic charges and Mandelstam kinematic variables, 
the decay distribution in terms of two lepton energies and three angles in the rest frame 
of the decaying neutralino is discussed. The consequences of the CP and CPT invariance 
on the polarized decay distributions are explained and new relations among the decay 
amplitudes, unique for Majorana particles, are derived. We illustrate in Sect. EJhow to 
probe the Majorana character and CP violation of the neutralino system through detailed 
analytical and numerical investigation of various observables in the rest frame of the 
decaying neutralino xH'- CP-even lepton energy /angular distributions, lepton invariant 
mass and opening angle distributions, and a CP-odd triple product of the neutralino spin 
vector and two lepton momenta. Finally, we summarize our findings and conclude in 
Sect. El 

*A clear independent evidence of the Majorana character of the neutralinos can be provided by an 
experimental identification of the selectron pair production in e~e~ collisions, which occur only via t— 
and it-channel neutralino exchange |1 51 117| . 
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2 Neutralino mixing 



In the MSSM, the four neutralinos Xi (i = 1> 2, 3, 4) are mixtures of the neutral U(l) and 
SU(2) gauginos, B and W 3 , and the SU(2) higgsinos, H® and H®. The neutralino mass 
matrix in the (B, W 3 , H±, H°) basis, 
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is built up by the fundamental SUSY parameters: the U(l) and SU(2) gaugino masses 
Mi and M 2 , the higgsino mass parameter fi, and the ratio tan/3 = v 2 /vi of the vacuum 
expectation values of the two neutral Higgs fields which break the electroweak symme- 
try. The existence of CP-violating phases in supersymmetric theories in general induces 
electric dipole moments (EDM). The current experimental bounds on the EDM's can be 
exploited to derive indirect limits on the parameter space |21j . 

By reparametrization of the fields, M 2 can be taken real and positive without loss 
of generality so that the two remaining non-trivial phases, which are reparametrization- 
invariant, may be attributed to Mi and [i: 



Mi = \Mi\e { 



(0 < < 2tt) 



(4) 



Since the matrix Ai is symmetric, one unitary matrix N is sufficient to rotate the gauge 
eigenstate basis (B, W 3 , H®, ) to the mass eigenstate basis of the Majorana fields Xi 



M 



diag 



N*MN* 



(5) 



The mass eigenvalues m ; (i = 1,2,3,4) in A^diag can be chosen real and positive by a 
suitable definition of the unitary matrix N. 

The most general 4x4 unitary matrix N can be parameterized by six angles and ten 
phases. It is convenient to factorize the matrix N into a diagonal Majorana-type M and 
a Dirac-type D component in the following way: 



N = MD 



(6) 



The matrix D can be written as a sequence of 6 independent two-dimensional rota- 
tions parameterized in terms of 6 angles and 6 phases. The diagonal matrix M = 
diag {e lQl , e"* 2 , e*" 3 , e iQ4 } is given in terms of Majorana phases a,i (0 < q.{ < it) |Hj. One 
overall Majorana phase is non-physical and, for example, ai may be chosen to vanish. 

Due to the Majorana nature of the neutralinos, all nine phases of the mixing matrix N 
are fixed by underlying SUSY parameters, and they cannot be removed by rephasing the 
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fields. CP is conserved if all the Dirac phases are mod it and the Majorana phases di 
mod 7r/2. [Majorana phases on = ±7r/2 describe different CP parities of the neutralino 
states.] As a consequence, all the matrix elements N ia are purely real or purely imaginary 
in the CP invariant case. 



Figure 1: Diagrams contributing to the leptonic three-body neutralino decay Xi ~^ 
X®£ + £~ ; the exchange of the neutral Higgs bosons is neglected because the contribution 
is strongly suppressed by the tiny electron and muon Yukawa couplings. 



3 Three— body leptonic neutralino decays 
3.1 Neutralino decay amplitude 

The diagrams contributing to the three-body leptonic decay process x® ~~ > are 
shown in Fig. ^ Here, the exchange of the neutral Higgs boson [replacing the Z boson] is 
neglected since the couplings to the first and second generation SM leptons, i = e and /x, 
are very small. In this case, all the components of the decay matrix elements are of the 
(vector-current) x (vector-current) form which, after a simple Fierz transformation of the 
slepton-exchange parts, may be written for the lepton final states as 



with the generalized bilinear charges D a p (a,/3 = L,R) for the decay amplitudes: 




d fc° - xV + r) 
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where the s-channel Z-boson, and the t- and u-channel slepton propagators are given 
by 

?7l- 

D 7 - 



s — m| + irrizT z 

D t L,R = ~ —2 



D U L,R — 2 r 1 r (^) 

in terms of the Mandelstam variables, s = (g_ + q + ) 2 , t = (pi — q~) 2 , and u = (pi — q+) 2 , 
where Pi, q + and g_ are the 4-momenta of the decaying neutralino x° and the positively 
and negatively charged leptons £ ± , respectively. The couplings Z^, g L ij and g Ri j are given 
in terms of the neutralino diagonalization matrix elements N ia (i, a = 1 - 4): 

% = \ (N l3 N; 3 - N ti N*,) , 

9Lij = 2 2 (N i2 c w + Nns w )(N* 2 c w + N^sw) 
4s w c w 

QRij = -^l^i (10) 

The complex couplings satisfy the hermiticity relations: 

Zij = Z*i, g L ij = glji, 9Rij = 9*Rji (11) 

so that, if the Z-boson and selectron widths are neglected in the Z and selectron prop- 
agators, all the bilinear charges D a p also satisfy similar relations. In the CP invariant 
case the mixing matrix elements Ni Q are purely real or purely imaginary implying that 
the couplings Z^, g L ij and g^j are also purely real or purely imaginary. However, these 
couplings are in general complex in the CP non-invariant case, having both non-trivial 
real and imaginary parts. 
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3.2 Neutralino decay distribution 

The absolute amplitude squared of the three-body leptonic decay x°i ~^ Xi^ + ^~ °f a 
neutralino Xi with its polarization vector n is given by [TH] (the full spin-density decay 
matrix can be found in Refs. |§1 122j) 

\V\ 2 (n) = 4 (m 2 - t) (t - m 2 )(Aq - N 3 ) + 4 (m 2 - u)(u - m 2 )(i\q + 7V 3 ) 
—8 mjWi! s + 16 mi (pinq_q + ) 
+8 (n • ?+) [m<(u - ^DOM + N! 3 ) - mi(m 2 - t)N^] 
+8 (n • g_) [rrii{t - m\) (N{ - N£) + m x {m 2 - u)JV£] (12) 

where n is the Xi s P m 4-vector and (pinq-q + ) = e^ upa p^n l 'q'Lq'^ with the convention 
^0123 = +1- For the sake of notation, we introduce the abbreviations, m, = m^o. The 
seven quartic charges A^ lt2 ,3,4 and iV{ 2 3 for the 3-body neutralino decays are defined in 
terms of the bilinear charges by 

Aq = ^ [|Dflfl| 2 + IDllI 2 + |Z^ L | 2 + |D L/? | 2 ] 

N a = i Ste (UiwUIfl + D LL D* RL ) 

N 3 = \ [\D LL \ 2 + \D RR \ 2 - \D RL \ 2 - \D LR \ 2 ] 

N A = ± 3m (D^D* B + D LL D RL ) 

N[ = \ [\D RR \ 2 + \D RL \ 2 - \D LR \ 2 - \D LL \ 2 ] 

N' 2 = l -Ke (D RR Dl R - D LL D* RL ) 

= 1 [\D RR \ 2 + \D LR \ 2 - \D RL \ 2 - \D LL \ 2 } (13) 

The quartic charges Aq are P-even, while the quartic charges Aq' 2 3 are P-odd. 

We choose the rest frame of the decaying neutralino as a reference frame to describe 
the 3-momenta of the decay products. The neutralino spin 3-vector ft = (0, 0, 1) defines 
the direction of the z-axis. Since the azimuthal angle is irrelevant, the 3- momentum vector 
of the negative lepton can be taken to fix the x-z plane and its polar angle is denoted 
by 9 (0 < 9 < 7r). The orientation of the neutralino decay plane (NDP) is then fully 
determined by specifying an additional angle a (0 < a < 2tt), so that by rotating the 
NDP by —a around the l~ direction it is brought to x-z plane, as depicted in Fig. |21 The 
differential decay distribution is written in terms of two dimensionless energy variables, 
X- = 2E e -/rrii and x + = 2E e +/rrii, and the two angles, 9 and a, as 

1 — T^~T — TT~ = ~uTT \ F o(x-,x+) + ($- ■n)F 1 (x-,x+) + (q+ ■ ft) F 2 (x_, x+) 
ax -dx+d cos 9 da lbir z L 
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+ n ■ (q_ x q + ) F 3 (x_,x+) 



(14) 



where cos 9 = q_ ■ n, q± = q±/\q±\ and lepton masses are neglected. The four kine- 
matic functions Fk(x-,x+) (k = 0-3) are expressed in terms of the dimensionless energy 
variables, x_ and x + , and the quartic charges as 

F (x_, x+) = (N t - N 3 ) + x + y + (Ni + N 3 ) - 2r fl (x_ + x + - 1 + rfj N 2 

Fi{x-, x+) = (N[ - N£ - r a x-x+ N' 2 

F 2 {x_,x + ) = -x + y + {N[ + N'^) + r tl x.x + N' 2 

F 3 (x-,x + ) = rnX-X + N A (15) 

where y± = 1 — x± — r\ x with rn = mx/rrii. All quartic charges are functions only of the 
energy variables x±, but independent of the orientation angles, 6 and a. 



x 



£+, x 




Figure 2: A schematic diagram of the kinematic configuration of the momenta and 
spin vector in the initial neutralino rest mass frame for the three-body neutralino decay 
Xi Xi^ + ^~ i n terms of two dimensionless energy variables x± and two angles 9 and a. 
The opening angle \ is uniquely determined by two variables x±. 

In the kinematic configuration of Fig. 121 with n = (0,0,1) taken along the positive 
z-axis and the negative lepton momentum on the x—z plane, the spin dependent scalar 
products are expressed in terms of the opening angle, Xi an d two angles, 9 and a, as 

g_ • n = cos 9, q + ■ h = cos \ cos ^ — srn X sm ^ cos a 

n ■ (q_ x q + ) = sin \ sin 9 sin a (16) 

where the cosine of two lepton momentum directions, cosx, is simply a function of two 
normalized energy variables, x_ and x + : 

cos x = 1 _ 2 (x- + x + — 1 + r^) / (x-X + ) (17) 
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and sinx = \/l — cos 2 \- The kinematically-allowed ranges for the angles, 6 and a, and 
the dimensionless energy variables, x + and are 

0<9, |<7T, 0<x ± <l-4, (l-x_)(l-x + )>r?!, x„+x+>l-4 (18) 

The allowed crescent-shaped x + and x_ region (see the left panel of Fig. [3J is restricted 
due to four-momentum conservation, while the full ranges of the angles 9 and a allowed. 

3.3 Majorana nature and implications of CP transformation 

Before presenting a few concrete numerical examples for probing the Majorana nature and 
CP violation in the neutralino system through the leptonic three-body neutralino decays, 
we investigate the implications of the invariance under CP and CPT transformations^ for 
the three body leptonic neutralino decays [HE]. 

In the rest frame of the decaying neutralino, three final particles form a decay plane and 
every three-momentum changes its sign under P transformation as well as T transforma- 
tion while C transformation exchanges the four momenta of two leptons. The polarization 
vector n does not change under P and C transformations, but it changes its sign under T 
transformation. Then, the CP operation transforms the momenta and spin vector in the 
decay processes as 

x± -> +x T , <f± -> -q T , q±-h^ -q T ■ n, n ■ (<?_ x q + ) -> -n ■ (g_ x q + ) (19) 
while the CPT operations transform the momentum and spin vectors as 

x± -> +x T , q± -> +q T , q±-n-^ -q T ■ n, n ■ (g_ x q + ) -> +n ■ (q_ x q + ) (20) 
with the condition that the complex conjugation of the decay amplitude is taken [23] . 

Since neutralinos are the Majorana particles, CP invariance in the three-body neu- 
tralino decay x° "~ > X^ + ^~ leads to additional relations among the bilinear charges defined 
in Eq.flSJ): 

d rl = mmD LL {t <-> u) (21) 

where rjij = ±i are the intrinsic CP parities [21] of the neutralinos, Xij> respectively, and, 
as a result, to the CP relations for the kinematic functions defined in Eq. (|15j) : 

F (x_,x+) = +F (x + ,x-) 

F 1 (x-,x+) = -F 2 (x + ,x_) 

Fs{x-,x + ) = -F 3 {x + , x.) (22) 

TThe naive time reversal transformation T reverses the direction of all 3-momcnta and spins, but 
does not exchange initial and final states. Quantities that are odd under CPT can be non-zero only for 
complex transition amplitudes with absorptive phases which can be generated, for example, by loops, 
and Breit-Wigner propagators. 
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in the three-body leptonic neutralino decays. We note that the CP relations (f2~2*|) can 
be satisfied only when the couplings Zy, guj and gmj are simultaneously purely real (for 
Vi = Vj = =t z ) or simultaneously purely imaginary (for rji = —r]j = ±i). 

On the other hand, CPT invariance leads to the relations among the bilinear charges: 

D LR = -D* RR (t «-> u) 

D RL = -D* LL {t <-> u) (23) 

These CPT relations are satisfied if the Z-boson and slepton widths are neglected, that is 
to say, if there are no absorptive parts in the process. In the approximation of neglecting 
particle widths, we have the following CPT relations for the kinematic functions: 

F (x_,x+) = +F (x + ,x_) 
Fi(x-,x+) = -F 2 (x + ,x-) 

F 3 (x_,x + ) = +F 3 (x + ,x_) (24) 

independently of the mixing character of neutralinos and whether CP is violated or not. 

CP-conserving absorptive parts appear in loop diagrams with on-shell propagators 
through final-state interactions or from the widths of intermediate unstable particles. 
In addition, CPT-odd asymmetries may arise from the interference between a dominant 
tree-level and a sub-leading loop diagram mediating the decay. However, those absorp- 
tive parts and interference effects are usually tiny in the leptonic decay involving only 
electroweak interactions, at most at a level of a couple of percents. In this light, we will 
ignore all the width effects and electroweak loop corrections in the following analytic and 
numerical analyses. 

4 Numerical Analyses 

In the numerical analyses below we adopt an mSUGRA scenario defined by 

m = 150 GeV, m 1/2 = 200 GeV, A = -650 GeV (25) 

at the GUT scale requiring the pole mass of the top quark m t = 178 GeV, and 

tan/3 = 10, sgn(/i)>0 (26) 

at the electroweak scale at which all parameters are derived with the RGE code SPheno 
[2*K] (very similar results are obtained with other RGE codes; for comparison of different 
codes see For the light neutralino and chargino masses we find 

m^o = 78.1 GeV, m^o = 148.5 GeV, m-± = 148.4 GeV (27) 
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and we note that the light neutralino and chargino masses and mixing angles are repro- 
duced fairly well with the tree-level formulae taking 

M 1 = 80GeV, M 2 = 158GeV, |/i|=415GeV, $ M = 0; tan/5 = 10 (28) 

The selectron and sneutrino masses derived with the RGE code SPheno are 

m e - L = 207.7 GeV, m e - fl = 173.1 GeV, m* e = 192.1 GeV (29) 

Relevant for our analysis, the derived branching ratios for leptonic three-body decays of 
the second lightest neutralino are 

Mxl - X?e + e") = 4.5%, Br(x° -> X>V) = 4.6% (30) 

and for the ei — > X2 e decay 

Br(g L -> x°e) = 28.4% (31) 

while other decay modes have Br(x° ~~ ► Xi r+,r ~) = 58.9%, Br(x° ~ ► XiQO) = 9.1%, 
Mxl -» $W) = 23.7%, Bi(e L -> x°e) = 21.4%, Br(£ L -> = 50.3%, Br(z> e -> 

X» = 19.6%, Br(z> e -> xf e) = 45%, and Brfc -> &xf) = 98.3%, Br(ti -> c^) = 1.6%. 

The production cross sections for selectron-pair production processes with unpolarized 
e + e~ beams at y/s = 500 GeV are as follows 

a{e+e^} = 273.4 fb, a{e%e^} = 113.5 fb, a{e+e^} = 80.7 fb (32) 

so that large ensembles of events, ~ 2 x 10 5 events for e^ef and at integrated 

luminosity of 1000 fb -1 , will be generated. Given the branching fractions in Eqs. ()30|) and 
fl3*Tj) . a sufficient number of events for the decays x° — > Xi e+e ~ an d X2 ~^ XiA* + A* - are 
expected to be selected*, allowing the analysis of the properties of the neutralino decay 
at great detail. 

For completeness, we note that other production channels for SUSY-related processes 
that are either open or with cross sections greater than 1 fb at y/s = 500 GeV have 
the following cross sections: cr{x\x\} = 264.8 fb, cr-fx^X^} = 159.1 fb, cr-fx^X^} = H6.8 
fb, <?{xlXi} = 294 fb, = 39.8 fb, a{/i+/i7j = 22.3 fb, a{f+f{} = 52.4 

fb, a{fff^} = 5.2 fb, a{f 2 + f 2 ~} = 17.1 fb, a{u e u*} = 662.2 fb, = 15.2 fb, 

<j{v T i>*} = 16.8 fb, cr{tit*} = 45 fb and a{h°Z} = 64.6 fb. Most of these processes can 
be separated by simple kinematical cuts. Note that the i> e v* e production process might 
be exploited for our purposes if the decay mode xf e °f m one hemisphere could be 
used to tag the decay mode y^u of the second electron-sneutrino in the other hemisphere. 
However, since the detailed experimental simulations to assess this possibility are beyond 
the scope of the paper, we do not include the u e u* channel in our analyses. 

*In Monte Carlo simulations wc conservatively assume that at least 1000 neutralino decay events can 
be selected and we evaluate all the relevant physical quantities with their tree-level formulas based on 
the parameter set <|28ll . 
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4.1 Lepton energy distribution 



The polarization-independent kinematic function Fq is symmetric with respect to the 
energy variables x + and x_ exactly in the CP invariant case and to a good approximation 
in the CP non-invariant case. One of the most decisive ways for confirming the Majorana 
nature of the neutralinos, namely the observation of the symmetric distribution of events 
on the (x-,x + ) Dalitz plane [I], can therefore be realized in the cascade decay process of 
Eq. (j2J) in which the x° res t frame is reconstructable. 

The left panel of Fig. El shows the Dalitz plot of the leptonic three-body decay X2 ~~ * 
Xi£ + £~ in the (x-,x + ) plane for the parameter set (}2*%)) with $1 = 0. The right panel 
shows the numbers of piled-up events with sign(x_ — x + ) = — (left histogram) and 
sign(x_ — x + ) = + (right histogram), simulated with 1000 events by a Monte Carlo 
method. We note that the difference of two numbers of events AiV ev m 24 is within the 
expected statistical error of AN cxp = y/N cv ~ 32 with N cv = 1000. 



Dalitz Plot 
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Figure 3: Left: the Dalitz plot of the neutralino decay X2 ~^ Xi^ + ^~ the (x-,x+) 
Dalitz plane. Right: the number of events with sign(x_ — x + ) = — (left histogram) and 
sign(x_ — x + ) = + (right histogram) . $1 = is taken and the parameter set i2S\) is used 
for the other relevant SUSY parameters. The event distribution is generated with the total 
number of events of 1000 by a Monte Carlo method. 



4.2 Lepton angular distribution 

Another kinematic distribution of great interest is the lepton angle distribution with 
respect to the neutralino polarization vector; the lepton angle distribution with respect to 
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the beam direction in the e + e~ reference frame has been discussed in jS]. Defining 0± to 
be the polar angle between the £ ± momentum and polarization vector n, the normalized 
lepton angle distribution can be written as 



1 dT 1 M 

1 acoso± 2 



(33) 



with cos^i = q± ■ h. In the approximation of all particle widths neglected, the slope 
parameters have to be equal, = £+, irrespective of whether the theory is CP invariant or 
not. This is the consequence of CPT invariance of the decay distribution of the Majorana 
particle, cf. the second relation of Eq. (|24|). As a result, the sum of two lepton angle 
distributions has to be flat, that is to say, independent of the polar angles. This is one of 
the genuine tests of the Majorana nature of the neutralinos. 
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Figure 4: Left: The normalized lepton angle distribution j j.V.Vj) of the neutralino decay 
X2 ~^ Xi^ + f-~ i the black solid line is for the negative charged lepton and the red dashed 
line for the positive charged lepton. Right: The $1 dependence of the slope parameter £_ 
for the parameter set 



The left panel of Fig. 0] shows the lepton angle distribution for the parameter set 
with the phase $1 = 0. The solid line is for the cosine of the negatively-charged 
lepton angle, cos#_, and the dashed line for the cosine of the positively-charged lepton 
angle, cos# + . A simple numerical analysis based on the number of events of iV ev = 
1000 shows that the CPT relation can be confirmed within 1— tr statistical uncertainty 
of about 8% for the whole range of cos6*±. It increases to about 10% for the range 
with I cos 6± I > 0.8 vetoed if the cut turns out necessary to avoid distortions of the 
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£ ± distributions differently by experimental selection criteria ^Hj- Certainly, it will be 
necessary to perform a more detailed and realistic experimental analysis including all the 
systematic uncertainties before drawing a more concrete conclusion. Such a comprehensive 
analysis is, however, beyond the scope of the present theoretical investigations. 

The quantities £±, denoting the slope of the lepton angle distribution with respect 
to cos#±, depend on the values of the SUSY parameters. In the right panel of Fig. 0] 
the $1 dependence of the CP-even quantities £± is shown demonstrating that the £± 
measurement can provide information on the CP-violating phase. 

4.3 Lepton invariant mass and opening— angle distribution 

The invariant mass mu of two final-state leptons in the decay x\ ~ * Xi^ + ^~ is a Lorentz- 
invariant kinematic variable so that it is straightforward to reconstruct the quantity ex- 
perimentally by measuring the four momenta of two final-state leptons in the laboratory 
frame. Furthermore, the distribution for the invariant mass mu is independent of the 
specific production process for the decaying neutralino, because the invariant mass distri- 
bution does not involve any spin correlations between production and decay [Sj. 

Near the maximum end point of the lepton invariant mass distribution, the neutralino 
Xi is produced nearly at rest, the Mandelstam variable s approaches to (m,2 — mi) 2 and 
the variables t and u become identical ~ mim 2 . Ignoring the particle widths in the 
propagators we find from Eq. (JZJ) that the decay amplitude can be written approximately 

as 

V ~ —u(Xi) \l [i*m(X L ) + 7 5 9fo(X L )] + R [iQm{X R ) + 75^e(X fl )] ) u(x° 2 ) (34) 

where Xl = Dll and Xr = D^r near the invariant mass end point and L^/R^ = 
, u(£~)7 /i PL ) ij v are the left-handed/right-handed lepton vector currents. The approx- 
imate form of the decay amplitude leads to the absolute amplitude squared of the following 
approximate form: 

\V\ 2 ~ r 21 (1 - r 21 ) 2 { MXl? + MXr) 2 } + 0{(3 2 ) (35) 

near the end point with the neutralino xl velocity j3 ~ a/1 — fin for the normalized 
invariant mass, fin = m/z/(m 2 — mi). Therefore, the invariant mass distribution exhibits 
a characteristic steep S"-wave (slow P-wave) decrease proportional to (3 (/3 3 ) when the CP 
parities of two neutralinos are the same (opposite), that is to say, if Xl,r is purely real 
(purely imaginary). On the other hand, in the CP non-invariant case, where both the real 
and imaginary parts of X^ r are non-vanishing, the invariant mass distribution decreases 
steeply in S'-wave 0. 

The threshold behavior of the invariant mass distribution can be understood by in- 
vestigating the selection rule of the orbital angular momentum by CP symmetry. In the 
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non-relativistic limit of two neutralinos, the orbital angular momentum L of the final 
two-lepton and LSP system satisfies the CP relation 



1 = -V1V2 (~1) L 



(36) 



With the same (opposite) relative CP parity leading to 771772 = —1 (+1), the selection rule 
(I36J1 forces the final two-lepton and LSP system to have L — (1) near the threshold. 
In other words, the neutralino axial-vector (vector) current corresponds to S*-wave (P- 
wave) excitations. 

Additional clear signature of the selection rule (|3l)j) is provided by the decay distri- 
bution with respect to the opening angle x between two leptons [Sj. Since the relation 
between the invariant mass and the opening angle is given by 



the invariant mass vein takes its maximum value for x = 71 f° r given x+ and x_, i.e. 
when the momentum directions of two leptons are opposite. In this case, the helicities 
of two leptons coupled to a vector current are opposite, rendering its total spin sum to 
be unity along the flight direction. Therefore, angular momentum conservation forces the 
orbital angular momentum to be zero. As a consequence, when two neutralinos have the 
same (opposite) CP parity, the opening angle distribution is enhanced (suppressed) near 
cosx = — 1- A similar argument based on total angular momentum conservation can be 
applied to show that the opening angle distribution exhibits the opposite property for 
cosx = 1, i-e. for the vanishing opening angle, x = 0. 

Based on the parameter set ()28|). we show in Fig. |S]the lepton normalized invariant 
mass distribution (left panel) and the opening-angle distribution (right panel) for $! = 
(black solid line) and $1 = it (red dashed line). The histograms are based on 1000 events 
generated by Monte Carlo. The invariant mass distribution decreases steeply (slowly) 
near the end point and the opening angle distribution is strongly enhanced (suppressed) 
near x = n f° r t wo neutralinos of the same (opposite) CP parities in the CP invariant 
case. Note however, that since in both cases the invariant mass distribution vanishes 
at the end point, the distinction between the S- and P-wave behavior might be tricky. 
On the other hand, the opening angle distribution can have a finite value (depending on 
the relative CP parity) facilitating the discrimination. In our numerical case, the x 2 / n 
(with n = 20 degrees of freedom) for the fits with $1 = O/71", respectively, are as follows: 
0.86/0.64 for the lepton invariant mass distributions (left panel), and 0.83/0.54 for the 
opening angle distributions (right panel), indicating that the theoretical curves fit the 
Monte-Carlo generated histograms with very good precision. Consequently, the invariant 
mass and/or opening angle distributions can provide us with a very powerful handle for 
determining the relative CP parity of two neutralino states xl an d 




(37) 
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Figure 5: The lepton normalized invariant mass distribution (left panel) and the opening 
angle distribution (right panel) in the three-body leptonic neutralino decay x® ~ * Xi^ + ^~ ■ 
The black solid lines are for $1 = 0, (i.e. neutralinos with the same CP parities) and the 
red dashed lines for $1 = it, (i.e. neutralinos with the opposite CP parities). 



4.4 CP— odd triple spin /momentum product 

The CP and CPT relations in Eqs. f!22l) and (|24jl enable us to construct a genuine CP-odd 
and CPT-even distribution^: 

F CP (x_,x+) = l - [F 3 (x_,x + ) +F 3 (x + ,x_)] (38) 

A typical CP-odd observable related to the CP-odd distribution (}3*Hj) is the triple product 

C P = n ■ (q + x q_) (39) 

formed with the spin vector h and two final lepton momentum directions. With the 
CP-odd observable we can construct the CP-odd asymmetry as follows: 

_ jV(0 C p > 0) - iV(Qcp < 0) = Iv i sinxF C p(x-,x+) dx.dx+ 
CP ~ N(0 CP > 0) + iV(0 C P < 0) ~ j v F (x_, x+) dx.dx+ [ ' 

where T> denotes the kinematically allowed (x_,x + ) Dalitz region defined in Eq. (|18|). 



^In the presence of absorptive CP-conserving phases due to particle widths or radiative corrections, 
there exist two additional CP-odd distributions; F (x-,x + ) — F n (x + ,x-) and Fi(x-,x + ) + F 2 (x + , 
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t\ 2 I ... i ... i ... i ... i ... i ... i ... i ... i ... i .. . 
0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Figure 6: The $i dependence of the CP-odd and CPT-even asymmetry Aqp of the triple 
scalar product for the parameter set Wty . We note that the 1-a statistical uncertainty of 
the CP-odd asymmetry is {I — A^ p ) /N cv ~ 3.1% with the number of events of N ev = 
1000. 

As described in Sect. CP violation in the neutralino system arises when the phases 
of Mi and/or /i are different from and/or it. These phases generally lead to large 
contributions to the EDMs. For the selectron masses under consideration in the present 
work, the experimental bounds on the electron EDM put the strongest constraints on the 
phases $i and $ M . Nevertheless, it has been demonstrated in a recent work [2^ that if 
the phase ^ is set to be or n, the full range of $i is allowed. In this light, we set $ M 
to be zero and vary $i in estimating the CP-odd asymmetry Aqp- 

Figure El shows the $i dependence of the CP-odd and CPT-even asymmetry Acp of 
the triple scalar product. Numerically we find that the magnitude of the asymmetry can 
be as large as 15% for the parameter set which is expected to enable us to clearly 
measure CP violation in the neutralino system directly with its expected 1-er statistical 
uncertainty of — Aq P )/N cv ~ 3.1% for the number of events of N ev = 1000. 

We note finally that if the measured value of Aqp turns to be very small, close to zero, 
the two-fold ambiguity $i = or $i = tt, can be resolved with the help of either the 
slope parameters £± or the lepton invariant mass/opening angle distribution, see Fig. 0] 
and Fig. respectively. 
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5 Conclusions 



Taking into account the possibility of having highly polarized neutralinos X2 an d of re- 
constructing their rest frames, we have performed a systematic analysis of the polarized 
neutralino decay X2 ~~ * Xi^ + ^~ i n its decay rest frame for all the physical implications due 
to the Majorana nature as well as CP violation of the neutralino system in the MSSM. 
We have demonstrated that the decay process {1} in the cascade channel like (J2J can be 
used to provide alternative powerful methods for probing the neutralino system in detail. 

In the CP invariant case, the Majorana nature of the neutralinos can be checked 
through the leptonic decay x® ~^ Xi^ + ^~ by verifying that in the x® res ^ frame: 

• the charged lepton energy distribution is identical irrespective of the electric charge 
of the lepton, 

• the sum of the negative and positive lepton angle distributions is independent of 
the lepton angles with respect to the neutralino polarization vector, 

• the relative CP parity of two neutralinos can be identified by measuring the threshold 
behavior of the invariant mass distribution near the kinematic end point and/or the 
dependence of the decay distribution on the opening angle near the angle close to 
7r and/or 0. 

In the CP non-invariant case, if all the absorptive parts (which are expected to be usu- 
ally tiny for such an electroweak decay) are ignored, both the lepton energy and angle 
distributions in the decay x[> — * Xi^ + ^~ can allow us to probe the Majorana nature of the 
neutralinos as in the CP invariant case. In addition, the three-body leptonic neutralino 
decays allow us to construct a CP-odd triple scalar product of the neutralino polarization 
vector and two lepton momenta. Numerically, we have found that for the parameter set 
f!28|) the CP asymmetry related to the triple product could be of the order of 15%, while 
satisfying the severe EDM constraints on the SUSY parameters. 

Finally, we emphasize that the analyses in the present work are only statistical and 
based on the assumption that both the degree of neutralino polarization and the efficiency 
of reconstructing the neutralino rest frame are 100%. Clearly, it is necessary to perform 
further detailed experimental analyses with realistic values of reconstruction efficiencies 
and with background processes included, which is beyond the scope of the present theoret- 
ical analyses. However, we think that the results of our theoretical studies are encouraging 
enough to motivate further detailed experimental investigations to assess the feasibility 
for probing the Majorana nature of neutralinos and CP violation in the neutralino system 
at future colliders. 
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